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Lattice studies of the infrared behavior of gluon and ghost propagators are a key probe of confine- 
ment scenarios in Yang-Mills theories. However, finite-volume effects become an important issue 
as the infrared limit is approached. By considering general quantities — namely an associated sus- 
ceptibility in the gluon case and properties of the lowest-lying eigenmode of the Faddeev-Popov 
matrix in the ghost case — one can derive rigorous upper and lower bounds for the propagators. 
The bounds allow a better control over the extrapolation of lattice results to the infinite-volume 
limit. In the case of the gluon propagator, an intuitive statistical argument suggests a precise vol- 
ume dependence for the bounds. This dependence is nicely confirmed by the lattice data, leading 
to a finite gluon propagator at zero momentum. At the same time, an enhancement of the ghost 
propagator in the infrared limit seems unlikely. Our analysis is applied to the case of Landau 
gauge and SU(2) gauge group, using the largest lattice sizes to date. 
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1. Introduction 

Confinement in QCD and in Yang-Mills theories in general is associated with long-range (in- 
frared) effects. It is thus necessary to study the infrared behavior of the theory's Green's functions 
using nonperturbative methods. These studies include numerical (lattice) as well as analytic meth- 
ods and consider basic (gauge-dependent) quantities — such as gluon and ghost propagators — 
in order to test the predictions of the so-called confinement scenarios (see e.g. [|I|] and references 
therein). In the case of lattice simulations, one has at one's disposal a true first-principles method, 
with no uncontrolled approximations. On the other hand, extreme care must be taken to extract the 
true infrared behavior of the propagators from lattice data, since significant systematic errors may 
affect the extrapolations that are needed in order to get physical results. The most important such 
errors are Gribov-copy effects 1 — related to the fact that the relevant objects are gauge-dependent 
quantities — and finite-size effects. The latter is especially important in the investigation of the 
infrared limit, since the smallest nonzero momentum that can be represented on a lattice of linear 
extension L is proportional to l/L. Thus, a sensible range of small momenta can only be properly 
simulated on a very large lattice. Here, we consider carefully the elimination of finite-size effects 
through a better control of the extrapolation of our data to the infinite-volume limit, as described 
below. 

The extrapolation of gluon- and ghost-propagator data to infinite lattice volume is a delicate 
task, since the correct volume dependence of the data may not be easily inferred from the behavior 
on medium-size lattices and since some quantities, such as the zero-momentum gluon propaga- 
tor, are quite noisy. For these reasons it proves very helpful to obtain constraints on the infrared 
behavior of the propagators, as the upper and lower bounds presented here. We remark that these 
bounds are valid at each lattice volume V and must be extrapolated to infinite volume, just as for the 
propagators. The additional advantage, besides establishing a range of allowed values for the prop- 
agators, is that the bounds are written in terms of "friendly" quantities — i.e. easier to compute, 
better behaved or more intuitive than the propagators themselves. It will therefore be more conve- 
nient to study the volume dependence of the bounds first, in order to assess the volume dependence 
of the propagators. We describe and apply the gluon and ghost bounds — in pure SU (2) theory 
and Landau gauge — respectively in Sections |2] and |3| below. As can be seen from our analysis, 
we obtain a finite nonzero gluon propagator and a tree-level-like ghost propagator in the infrared 
limit. Possible implications of these results for the currently accepted confinement scenarios are 
discussed in the Conclusions. 



2. Gluon Bounds 

Rigorous upper and lower bounds have been inuoduced in Ref. [|]] for the gluon propagator at 
zero momentum, defined as 

D(0) = 5(«FT)M <0)|2> - <21) 



1 The problem of Gribov-copy effects has been extensively studied on the lattice |g||3| ^, g] . We comment briefly on 
this issue in our Conclusions. 
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where A b „{p) is the Fourier transform of the gluon field A b „(x) in pure SU(N C ) gauge theory, (•) 
stands for the path integral (Monte Carlo) average, V = N d is the lattice volume and we consider d 
space-time dimensions. Let us define the quantity 

M(0) = 5(]vFT)£l A> >l- i22> 

It is straightforward to show that this quantity is related to D(0) as 

V (M(0)) 2 < D(0) < Vd(N 2 - 1) (M(0) 2 ) , (2.3) 



which provides us with rigorous upper and lower bounds for D(0) that must be satisfied at every 
volume V. We can now try to interpret the quantities (M(0)} 2 and (M(0) 2 ), to obtain perhaps an 
understanding of their volume dependence. We start by noting that if we take the above "magneti- 
zation" without the absolute value, i.e. considering 

M ' (0> = TwhuP'^ ■ (24) 



we get a null Monte Carlo average: (M'(0)) = 0. Because of the absolute value, the quantity defined 
in Eq. ( ^2| ) has a nonzero average at finite V, but it should go to zero at least as fast as V~ x l d ', as 
shown in [{7|]. We now note that V{M(0) 2 ) is essentially the susceptibility associated with the 
magnetization M'(0) (since the average of this magnetization is zero). For a (i-dimensional spin 
system one thus expects to see V (M(0) 2 } ~ const, i.e. the statistical variance of the magnetization 
is proportional to the inverse of the volume, a behavior known as self-averaging. At the same time, 
considering the statistical fluctuations in the Monte Carlo sampling of M(0), we would expect 
(M(0)) 2 to have the same volume dependence as (M(0) 2 ) [0]. 

The simple statistical argument presented above suggests that both (M(0)} 2 and (M(0) 2 } 
should show a volume dependence as 1/V, implying (for d > 2) a much stronger approach to 
zero than the limiting behavior for M(0) obtained in [[7|] and mentioned earlier. On the other hand, 



the suppression with 1 /V is compensated by the volume factor for both bounds in Eq. (23). Con- 
sequently, if this suggested behavior for the susceptibilities is verified, D(0) converges to a nonzero 
constant in the infinite-volume limit. Note that the bounds in Eq. ( pT3| ) apply to any gauge and that 
they can be immediately extended to the case D{p) with p ^ 0. 

We have investigated the volume dependence of the bounds for pure SU (2) gauge theory in 
Landau gauge, considering physical lattice volumes of up to a V « (27 fm) . We find remarkably 
good agreement with the predicted 1 /V behavior for (M(0)) 2 and (M(0) 2 ), as can be seen in plots 
and tables in Ref. More precisely, by fitting the two quantities to \/V a we get the exponents 
a respectively 0.995(10) and 0.998(10). Analogously, an analysis for the SU (3) case (considering 
somewhat smaller volumes) yields the exponents 1.058(6) and 1.056(6) A similar behavior 
is also obtained by a study introducing a modified gauge-fixing procedure (in order to check for 
possible Gribov-copy effects) [0]. Finally, a finite nonzero gluon propagator has been recently 



obtained using improved actions and anisotropic lattices [10]. We remark that this behavior has 



also been clearly observed on very large lattices in 3d [11, 6|] but not in 2d [12, §]. 
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3. Ghost Bounds 



Rigorous lower and upper bounds for the ghost propagator G(p) were proposed in [ 13]. We 
recall that G(p) is given by the inverse of the Faddeev-Popov (FP) matrix ^# and that an infrared 
enhancement of G(p) with respect to the tree-level ghost propagator G(p) ~ p~ 2 is generally ex- 
pected as a sign of confinement. By straightforward calculations — independent of the ones per- 
formed in the gluon case — we can establish bounds for the ghost propagator. In Landau gauge, 
for any nonzero momentum p, one finds 

JJYZTT T~ I- \ Wmin{a,p)\ 2 < G(p) < ^- , (3.1) 

jv c 1 "mm a "mm 

where A min is the smallest nonzero eigenvalue of the FP operator and y m i,,(a,p) is the corre- 
sponding eigenvector. Note that the upper bound is independent of the momentum p. If we now 
assume A m! „ ~ L~ v and G{p) ~ p- 2 - 2K a t small p, we have that 2 + 2k < v, i.e. v > 2, is a neces- 



sary condition for the infrared enhancement of G(p). A similar analysis can be carried out [ ]14[ ] for 
a generic gauge condition. Consider the Gribov region £2, where all eigenvalues of ^# are positive. 
In the infinite- volume limit, entropy favors configurations near the Gribov horizon d£l (where A m ;„ 



goes to zero). Thus, inequalities such as (3. 1 ) can tell us if one should expect an enhancement G{p) 



when the Boltzmann weight gets concentrated 2 on dQ.. 



Our study in the SU (2) Landau case [13], using the very large lattices mentioned in the pre- 
vious section, suggests that v < 2 (for d = 4). This tree-level-like behavior is confirmed if one 
considers the dressing function p 2 G(p). Indeed, the data for the dressing function can be well 
fitted by a — Mog(l +cp 2 ) [|13||, supporting fc = 0. This is also observed in d = 3. For d = 2 
enhancement is observed, with a behavior ~ p~ lK and K* between 0.1 and 0.2 [12, [Oy . 



4. Conclusions 

By using rigorous upper and lower bounds to constrain the infrared behavior of gluon and ghost 
propagators, we obtain a finite nonzero gluon propagator at zero momentum and an essentially con- 
stant ghost dressing function in the infrared limit. These results seem to contradict the commonly 



accepted confinement scenarios of Gribov-Zwanziger and Kugo-Ojima [16]. However, as pointed 



out in [17], the above results are not completely in disagreement with the Gribov-Zwanziger ap- 
proach. In particular, it has been recently shown JT8| ] that using the Gribov-Zwanziger approach, 
i.e. by restricting the functional integration to the Gribov region Q., one can also obtain in 3d and 
4d a finite nonzero gluon propagator and a tree-level-like ghost propagator in the infrared limit. It 



is interesting that the same approach cannot be applied to the 2d case Q19p. Let us also note [|1J 
that even though the Gribov-Zwanziger and the Kugo-Ojima confinement scenarios seem to predict 
similar infrared behavior for the propagators, it is not clear how to relate the (Euclidean) cutoff at 
the Gribov horizon to the (Minkowskian) approach of Kugo-Ojima [p0|]. 

Similar results for the gluon and ghost propagators are obtained by various groups using very 



large lattice volumes [21], both in the SU (2) and in the SU (3) cases. [The equivalence between the 



For example, in Ad Maximally Abelian gauge one sees that goes to zero at large volume but the ghost propa- 
gator stays finite at zero momentum |p^]. 
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infrared propagators in SU (2) and SU (3) gauge theories can be seen e.g. in [22].] Of course, one 
should also recall that the region Q. is actually not free of Gribov copies and that the configuration 
space should be identified with the so-called fundamental modular region A. On the other hand, 
the restriction to A and the numerical verification of the Gribov-Zwanziger scenario are separate 
issues [|l7|]. Indeed, this scenario is based on the restriction of the configuration space to the region 
£2, which includes A. Finally, as explained in the restriction to A can only make the ghost 
propagator less singular, as confirmed by recent lattice data [^J. 
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